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ABSTRACT 

IRAS 09104+4109 is a rare example of a dust enshrouded type 2 QSO in the centre of a 
cool-core galaxy cluster. Previous observations of this z=0.44 system showed that as well 
as powering the hyper-luminous infrared emission of the cluster-central galaxy, the QSO is 
associated with a double-lobed radio source. However, the steep radio spectral index and mis- 
alignment between the jets and ionised optical emission suggested that the orientation of the 
QSO had recently changed. We use a combination of new, multi-band Giant Metrewave Radio 
Telescope observations and archival radio data to confirm that the jets are no longer powered 
by the QSO, and estimate their age to be 120-160 Myr. This is in agreement with the ~70- 
200 Myr age previously estimated for star-formation in the galaxy. Previously unpublished 
Very Long Baseline Array data reveal a 200 pc scale double radio source in the galaxy core 
which is more closely aligned with the current QSO axis and may represent a more recent 
period of jet activity. These results suggest that the realignment of the QSO, the cessation of 
jet activity, and the onset of rapid star-formation may have been caused by a gas-rich galaxy 
merger. X-ray observations reveal a spiral structure in the ICM which suggests that the cluster 
is in the process of relaxation after a tidal encounter or merger with another system; such a 
merger could provide a mechanism for transporting a gas-rich galaxy into the cluster core 
without stripping its cold gas. A Chandra X-ray observation confirms the presence of cavi- 
ties associated with the radio jets, and we estimate the energy required to inflate them to be 
^7.7 x 10 60 erg. The mechanical power of the jets is sufficient to balance radiative cooling in 
the cluster, provided they are efficiently coupled to the intra-cluster medium (ICM). We find 
no evidence of direct radiative heating and conclude that the QSO either lacks the radiative 
luminosity to heat the ICM, or that it requires longer than 100-200 Myr to significantly impact 
its environment. 

Key words: galaxies: clusters: individual (CL 09104+4109) — galaxies: individual 
(IRAS 09104+4109) — galaxies: quasars: general — X-rays: galaxies: clusters — galaxies: 
active — cooling flows 



1 INTRODUCTION 

It is now widely agreed that feedback from active galactic nuclei 
* E-mail: ejos@star.sr.bham.ac.uk (AGN) is likely to be the dominant process governing star forma- 
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tion and gas cooling in systems ranging from large galaxies to the 
most massive galaxy clusters ( |McNamara & Nul sen 2007 ; [Peter-| 
|son & Fabian 2006). In the nearby Universe there is a wealth of 
evidence that radio galaxies in groups and clusters have an im- 
pact on the surrounding intra-cluster medium (ICM) by inflating 
caviti es (e.g.,|Rafferty et al.|2006[ |Birzan et al.|2008[ |Dunn et al.| 
2010| [O'Sullivan et al.||2011|, driving shocks and sound waves 
(e.g., |Nulsen et al.||2005| |Fabian et al.||2006| |Sanders & Fabian 
2007 IForman et al.|[2007>, and causing gas mi xing (e.g., |Fabian 



etaLl2 005 ; Wise et al. 2007 , Bla nton et alT2 009 ; Simionesc u et al. 



20091 |Kirkpatrick et al.|20TT) . 

While a detailed, universally accepted description of the 
mechanisms underlying the origin and power of AGN jets remains 
elusive, the combination of the results from numerous analytic 
studies and increasingly sophisticated numerical simulation stud- 
ies (e.g.,|Koide|2003||Gammie et al.|2004l|De Villiers et al.|2005| 



Hawley & Krolik 2006; Komissarov 2007; Punsly 2007 , Beckwith 
et al. 2008; McKinney & Blandford 2009) do suggest an emerg- 
ing consensus on a few key elements of the picture (c.f. |Benson &| 
|Babul|2009| for a more detailed overview): 

First, jets are an ubiquitous feature of magnetized accretion 
flows onto black holes (|De Villiers et al.|2005[[Punsly|2007[pcK^ 



inney & Blan dford 2009]). The jets may be highly relativistic, nar- 
rowly collimated and Poynting-flux dominated or broad, mildly rel- 
ativistic and matter-dominated. 

Second, the launching of powerful jets is most efficient when 
the accretion flow in the vicinity of the black hole has a large 
scale-height, as is the case when the flow is radiatively inefficient 
and advection-dominated (e.g., Meier 2001 ; Churazov et al. 2005 
|Nemmen et al.|2007] ). Radiatively inefficient flows are believed to 
arise at low accretion rates (i.e M < O.OlMEdd) (e.g., |Narayan &| 
|Yi|1994| >. Above this threshold accretion rate, the flow transitions 
to a radiatively efficient thin disk. There is a considerable body of 
work suggesting that the jet production is suppressed in thin disks 
( |Livio et al.|1999||Meier|2001|[Maccarone et al.|2003| l. At very high 
accretion rates (i.e M > 0.1 — 0.5MEdd), the nature of the flow 
and even its radiative properties are not well understood, though 
there are suggestions that as the flow becomes optically thick, the 
disk scale height increases and powerful jets may be launched. 

Third, the outflows draw their energy from both the gravita- 
tional energy released by the accretion flow as well as the rotational 
energy of the black hole itself. This link between black hole spin 
and outflow power has long been indicated from theoretical consid- 
erations (c.f .,|Blandford & Znajekp9T7l|Punsly & Coroniti|1990| 
|Meier| 1999 2001 it and in recent years, has received further support 
from numerical studies (c.f., Hawley & Krolik 2006 , Komi ssarov] 
|2007l |Punsly|[2007l |McKinney & Blandford||2009} . MHD simu- 
lations show that the magnetic fields anchored in the rotating ac- 
cretion flow will give rise to both an outflow of material from the 
accretion flow and the surrounding corona, as well as a hydromag- 
netic jet that originates in the innermost regions of the accretion 
flow. Of these two components, the former is thought to dominate 
if the black hole is spinning slowly; the frame-dragging of the ac- 
cretion flow within the ergosphere can greatly magnify the ampli- 
tude of the outflow ( |Meier|1999|[2Wl|Ne7nmen et al.|2007|[BerTl 
|son & Babul 2009). The power of the electromagnetic component 
rises very steeply at high black hole spins. When the accretion flow 
reaches the black hole event horizon, the gas is expected to drain 
off the field lines while the field lines, which for all intents and 
purposes can be thought of as being anchored on the black hole's 
event horizon, establish a structure similar to that first described by 
Blandford & Znajek ( 1977 ). In the event that the black hole is spin- 



ning, the winding of the magnetic field lines in the ergosphere will 
drive helical twists that manifests as highly collimated, Poynting- 
flux dominated jets. In the picture outlined here, jets arise regard- 
less of whether the black hole is spinning or not; however, a rapidly 
spinning black hole can greatly enhance the total outflow power. 

However, it has been argued that the mechanical power of ra- 
dio jets is unlikely to represent more than ~10% of the total energy 
budget of all types of AGN (Cattaneo & Best 2009 ). The remaining 
90% is released radiatively, generally by AGN with little or no jet 
activity, but it is as yet unclear whether this radiation can effectively 
heat gas on galactic or cluster scales. 

Current structure formation models suggest that feedback 
from black holes at redshifts > 1 is required to establish the shal- 
low gas entropy profiles in groups (e.g., McC arthy et al.||2011[ > 
and quench star formation in giant galaxies, so as to produce the 
distribution of luminosities and colours observed today ( |Bower| 
et al. 2006, Croton et al. 2006). It is widely assumed that radiative 
"quasar-mode" feedback, which should dominate at this epoch, 
can effectively heat the ICM via a wind which blows any cool gas 
and dust out of the galaxy (e.g., Ho pkins et al.| i2005 ) and drives 
shocks into the ICM. However this is a largely untested hypothe- 
sis, since powerful quasars are rare at low-redshifts (z < 1) where 
current instrumentation could resolve such winds, particularly in 
the dense cluster environments where their impact could be most 
easily studied. The handful of examples which have been exam- 
ined using Chandra and X MM-Newton (e.g., |Belsole et al.|2007| 
Siemiginowska et al. 2010, Russell et al. 2010, 2012) show no evi- 
dence of strong quasar-mode heating, but the difficulty of identify- 
ing disturbed structures or ascertaining the current accretion state 
means that the issue is as yet unresolved. 

We have chosen to investigate this issue by examining 
IRAS 09104+4109 (also known as IRAS J09 13454+405628) a 
quasar and FR-I radio galaxy located in the dominant galaxy of 
the 2=0.442 galaxy cluster CL 09104+4109 (also known as MACS 
J0913. 7+4056 and hereafter referred to as CL09). In this paper we 
combine new multi-frequency Giant Metrewave Radio Telescope 
(GMRT) and archival Very Large Array ( VLA) and Very Long Base- 
line Array (VLBA) observations with a reanalysis of the Chandra 
and XMM-Newton data, to examine in detail the interaction be- 
tween the AGN, radio jets and ICM in CL09. Section fTTTj describes 
our current knowledge of the cluster and AGN, based on previous 
multiwavelength studies in the literature. In Sections|2]and[3]we de- 
scribe the radio and X-ray observations, and their reduction. The 
resulting radio images are shown in Section [4] and X-ray imaging 
and spectroscopic analysis is described in Section [5] Section [^de- 
scribes the spectral analysis of the radio source and presents limits 
on its age and physical properties. In Section [7] we summarise our 
observational results, examine the impact of the radio source on 
the ICM, and consider the origin of various structures in the clus- 
ter. In Section [5] we discuss the likely history of the cluster, BCG 
and AGN, and the ability of the QSO to balance radiative cooling 
in the ICM. Our conclusions are summarised in Section|9] 

Throughout the paper we assume a flat cosmology with Ho = 
73 km s" 1 Mpc -1 , Q x = 0.73 and fi = 0.27. For the redshift of 
CL09 (z = 0.4418, Sloan Digital Sky Survey Data Release 3), this 
gives a luminosity distance 2372 Mpc, an angular diameter dis- 
tance 1140 Mpc, and angular scale of l"=5.5 kpc. We assume a 
galactic hydrogen column of 1.8 x 10 20 cm -2 , taken from the Lei- 
den/Argentine/Bonn Galactic Hi Survey ( Kalberla et al. 2005 ). 
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1.1 Summary of previous studies 



IRAS 09104+4109 is a hyper-luminous infrared galaxy (HLIRG, 
|Sanders & Mir abel 199 6|>, having an infrar ed luminosity 
>10 4B erg s" 1 {>W 6 L e , |Kleinmann et al.|l988| . The flux den- 
sity of the galaxy peaks at ~50 /im and spectral energy distri- 
bution (SED) modelling suggests that ~70% of the emission be- 
tween UV and Sub-mm is produced by a dust-shrouded active 
nucleus ( |Ruiz et aL]|2010) with a total bolome tric luminosity 



(2.3 - 3.8) x 10 4Y erg s 1 jVignali et al.' 



20111. The AGN is thus 



classed as a type II QSO, with dust reprocessing almost all its op- 
tical and ultraviolet emission into the infrared (Kleinmann et al.| 
[l9"88l|Vignali et al.|2011||Ruiz et al.|20Tol >. 

The remaining ~30% of the emission from the galaxy arises 
from star-formation ( |Vignali et aLl|2011| |Ruiz et al.||2010) . The 
galaxy colour profile becomes significantly bluer within 20 kpc of 
the centre, and the star-formation rate ( SFR) based on the e quiva- 
lent width of [Oil] is 41±12 M Q yr" 1 ( |Bildfell et al.|2008| . Mod- 



elling of optical and near-ultraviolet (NUV) photometry suggests 
that recent star formation with a probable age of ~70-200 Myr has 
contributed a small but significant fraction of the stellar mass in the 
brightest cluster galaxy (BCG, probably >5% |Pipino et al.[2 009). 

Although the QSO i s obscured, the BCG contains relatively 
little dust, <5x 10 7 M i Combes et al.|2011 1, with a mean tem- 



perature ~120 K i Kleinmann et al.| 19881. Based on a detection 
of CO emission Com bes et al.| ( [201*1 1 estimate that the BCG con- 
tains only ~3.2x 10 9 Mq of molecular hydrogen, significantly less 
than many star-formation dominated ULIRGs (see also Evans et al. 
1998). It is possible that the star-formation and QSO have con- 
sumed or evaporated much of the cool material in the system, or 
that the winds from the QSO have swept the BCG clean. 

Ground-based spectrophotometry shows a "plume" of [Om] 
emission extending ~5" (~27 kpc) north of the galaxy, wit h a ve 
locity only ~100 km s _1 different from that of the nucleus (jCraw 



|ford & Vanderriest 1996). Hubble Space Telescope (HST) imag- 
ing shows this to be an extended cool gaseous filament ( |Armus| 
|et al.| 1999| labelled "inner filament " in Fig. [7} illuminated by spec - 
trally hard photons from the AGN (Crawford & Vanderriest 1996), 
while ground-based spectroscopy suggests that it contains dust 
l |Tran et aL"1|2000| >. Other filamentary structures around the BCG 
include apparent gaseous filaments north and south of the BCG at 
radii of 5-9" (27-50 kpc, labelled "outer filaments") and smaller- 
scale "whiskers" extending 5-7 kpc from the BCG on roughly east- 
west axes l |Armus et al.||199^| see Fig. [TJ. These filaments have 
significantly bluer colours than the [Om] plume, but their origin is 
unclear. Possibilities include material stripped from cluster galax- 
ies, gas cooling from the hot ICM, or material illuminated by the 
AGN through openings in the dust shrouding the BCG. 

On arcsecond scales HST imaging shows highly polarized 
emission north and south of the AGN, with a polarization angle 
roughly perpendicular to its extension. These structures are thought 
to be the optical/UV scattering cones of the AGN iHin es et al.| 
1999). Ground-based polarimetry also shows the [Om] filament 
to be highly polarized ( |Tran et al. 2000). Attempts to determine 
the true opening angle of the cones and their inclination to the line 
of sight result in a range of solutions (half-opening angles 15°- 
40° and inclinations 34° -50° |Tran et al.|2000| |Hines et al.|1999| >, 
but the projected half-opening angle measured from the HST im- 
ages is 35°±10° pines et al.|1999l >. 

|Hines & Wilis] ( [1993) , using VLA 1.4 and 5 GHz observa- 
tions, found that IRAS 09104+4109 also hosts a double-lobed ra- 
dio source (NVSS 1091345+405630) with straight jets extending 



Figure 1. False colour HST PC2 image of IRAS 09104+4109 and its sur- 
roundings. An 814W filter image is used for the red band, a 622W image 
for green, and a scaled 622W subtracted from the 814W image for the blue 
band. This combination is chosen to emphasize the small-scale filamentary 
structures referred to in the text. Correcting for the redshift of CL09, the 
622W and 814W filters have central wavelengths (and widths) of ~4300A 
(~635A) and ~5550A (~1070A) respectively. The redder colour of the 
inner filament corresponds to strong [Om] 5007A emission from this struc- 
ture (Armus et al.| 1999) . 



~60 kpc northwest and southeast of the nucleus, whose radio lu- 
minosity places it in the FR I- FR II transition region. Its spectral 
inde>Qis 01.4-5 GHz ~ 1-4, fairly steep for a radio galaxy . The 
jet axis lies just outside the ionisation cones of the AGN, and it has 
been suggested that this misalignment, along with the brightness 
of the radio core (unresolved in 6cm and 20cm VLA observations) 
may indicate that AGN has recently changed orientation, with a 
new jet aligned close to the line of sight producing the bright core 
via relativistic beaming ( |Hines & Wil ls 1993; Hines et al.] |1999| 
|Armusetal.|1999| >. 

ROSAT HRI observations of the CL09 cluster showed ex- 
tended X-ray emission centred on the BCG, with a central hole 
~4" (~20 kpc) across <|Fabian & Crawford|l995). ROSAT imaging 



and ASCA spectral analysis both indicate the presence of a strong 
cooling flow in the cluster core ( |Fabian & Crawford|1995||Craw-| 
|ford & Vand erriest 1996). More recent Chandra observations have 
confirmed the presence of a strong cool core with a short cool- 
ing time (2.1 Gyr within ~50 kpc) and shown that the hole is one 
of two cavities located northwest and southeast of the core along 
the jet axis I Hlavacek- Larrondo et al.|2011||Iwasawa et al.|2001| >. 
The enthalpy of these cavities is more than sufficient to balance ra- 
diative losses within the cooling region I Hlavacek-Larrond o et al.| 
[20TT| hereafter HL11). 



1 The radio spectral index a is defined according to S v oc v 
is the flux density at the frequency v. 



1 , where S v 
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2 RADIO OBSERVATIONS AND DATA REDUCTION 

2.1 GMRT Observations 

CL09 was observed using the GMRT at 240 MHz, 330 MHz, 610 
MHz (project 15SRC01) and 1.28 GHz (project 17_050) in a full- 
synthesis run of approximately 11 hours at 1.28 GHz and 9 hours 
at each of the other frequencies (including calibration overheads). 
Details on these observations are summarised in Table ^ which 
reports: observing date, frequency and total bandwidth, total time 
on source, full-width half maximum (FWHM) and position angle 
(PA) of the full array, and rms level (la) at full resolution. 

The 330 MHz and 1.28 GHz data were recorded using the up- 
per and lower side bands (USB and LSB, respectively) simultane- 
ously, which provide a total observing bandwidth of 32 MHz. The 
240 MHz and 610 MHz observations were carried out in simul- 
taneous mode with 32 MHz bandwidth (USB+LSB) at 610 MHz 
and 8 MHz at 240 MHz. The data were collected using the default 
spectral-line mode with 128 channels for each band, resulting in a 
spectral resolution of 125 kHz per channel. The data sets were cali- 
brated and reduced using the NRACj^Astronomical Image Process- 
ing System (AIPS) package. The data were initially inspected using 
the task SPFLG to identify and remove bad channels and visibili- 
ties affected by radio frequency interference (RFI). The data were 
then calibrated. The flux density scale was set using amplitude cal- 
ibrators (3C48, 3C147 and 3C286), observed at the beginning and 
at the end of the observing run, and the Perley & Taylor ( 1999]> ex- 
tension of the (Baars et al. 1977) scale. The source 0834+555 was 
used as phase calibrator at all frequencies. 

The bandpass calibration was carried out using the flux density 
calibrators. A central channel free of RFI was used to normalise the 
bandpass for each antenna. After bandpass calibration, the central 
84 (48) channels were averaged to 6 channels of -~2 MHz (1 MHz) 
at 330 MHz, 610 MHz and 1.28 GHz (240 MHz) to reduce the size 
of the data set and, at the same time, to minimise the bandwidth 
smearing effects within the primary beam of the GMRT antenna. 

After further careful editing of the in the averaged data, a num- 
ber of phase-only self-calibration cycles and imaging were carried 
out for each data set. The large field of view of the GMRT required 
the implementation of wide-field imaging technique in each step of 
the data reduction, to account for the non-planar nature of the sky. 
The USB and LSB were calibrated separately. The final data sets 
were further averaged from 6 channels to 1 single channej^] and 
then combined together to produce the final images. These were 
corrected for primary beam pattern of the GMRT antenna using 
PBCOR in AIPS. The rms noise level (lcr) achieved in the final full 
resolution images are summarised in Table[T| 

The average residual amplitude errors on each individual an- 
tenna are < 5% at 610 MHz and 1.28 GHz and < 8% at 240 MHz 
and 330 MHz (e.g., |Chandra et al.|2004) . Therefore, we can conser- 
vatively assume that the absolute flux density calibration is within 
5% at 610 MHz and 1.28 GHz and 8% at lower frequency. 

2.2 VLA Archive Data 

We extracted all observations of CL09 useful for the analysis pre- 
sented in this paper from the VLA public archive and re-analysed 

2 National Radio Astronomy Observatory. 

3 Bandwidth smearing is relevant only at the outskirts of the wide field, 
and does not significantly affect the region presented and analysed in this 
paper. 



them. The observing details are provided in Table [2] which shows 
the project code and array configuration in the first two columns; 
the other columns provide the same information as Table [T] 

Calibration and imaging were carried out using AIPS, follow- 
ing the standard procedure, i.e., Fourier-Transform, Clean and Re- 
store. We applied phase-only self-calibration to remove residual 
phase variations and improve the quality of the images. Correction 
for the primary beam attenuation was applied to the final images 
using the task PBCOR in AIPS. All flux densities are on the Perley] 
|& Tayl or ( 1999 1 extension of the Baar s et al.H1977) scale. Average 
residual amplitude errors are < 5% at all frequencies. 



2.3 VLBA Archive Data 

CL09 was observed at 1.4 GHz by the VLBA on 2003 March 13 
(project BH01 10) for approximately 12 hours in phase referencing, 
with a recording bandwidth of 8 MHz. We extracted the calibrated 
data from the VLBA archive, processed with the VLBA data cali- 
bration pipeline (Sjouwerman et al. 20051 implemented in AIPS. 
The data set was further self-calibrated in phase with the AIPS task 
CALIB to reduce phase fluctuations, and used to produce the final 
image with a resolution of 8.3 mas x 6.3 mas. The rms noise level 
measured on the image plane is 25 pjy beam -1 . 



3 X-RAY OBSERVATIONS AND DATA REDUCTION 

CL09 was observed during Chandra cycle 10, on 2009 January 06 
for ~77 ks (ObsID 10445), with the ACIS-I instrument operating 
in very faint telemetry mode. A summary of the Chandra mission 
and instrumentation can be found in Weisskopf et al. (20021. The 
data were reduced and analysed using CIAO 4.3 and CALDB 4.4.3 
follow ing techniques similar to those described in O'Sullivan et al. 
1 2007} and the Chandra analysis thread^] The observation did not 
suffer from significant background flaring, and the final cleaned 
exposure time was 73.3 ks. 

Point sources were identified using the WAVDETECT task, with 
a detection threshold of 2.38 x 10 -7 , chosen to ensure that the task 
detects ^ 1 false source in the ACIS-I field, working from a 0.3-7.0 
keV image and exposure map. All point sources were excluded, ex- 
cept in some cases the source corresponding to the cluster-central 
AGN. Spectra were extracted using the SPECEXTRACT task. Spec- 
tral fitting was performed in XSPEC 12.7.0e. Abundances were 
measu red relative to the abundance ratios of Grevess e & Sauvall 
( |1998) >. A galactic hydrogen column of 0.018xl0 22 cm" 2 and a 
redshift of 0.44 were assumed in all fits. Spectra were grouped to 
20 counts per bin, and counts at energies outside the range 0.5-7.0 
keV were generally ignored during fitting. 

Background spectra were drawn from the standard set of CTI- 
corrected ACIS blank sky background events files in the Chandra 
CALDB. The exposure time of each background events file was 
altered to produce the same 9.5-12.0 keV count rate as that in the 
target observation. Very faint mode background screening was ap- 
plied to both source and background data sets. 

XMM-Newton has also observed CL09, for ~14 ks (Obsld 
0147671001). This observation is described in detail in |Piconcelli| 
|et al.| |2007) and a detailed summary of the XMM-Newton mission 



4 http://asc.harvard.edu/ciao/threads/index.html 
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Table 1. Details of the GMRT radio observations 



Observation 


Frequency 


Bandwidth 


Integration 


FWHM, PA 


rms 


date 


(MHz) 


(MHz) 


time (min) 


(full array, "x", °) 


(mJy beam -1 ) 


Dec 2008 


240° 


8 


390 


19.9 x 19.3, -7 


1.3 


Mar 2009 


327 


32 


280 


11.1 x 9.1, 17 


0.20 


Dec 2008 


610 a 


32 


390 


6.1 x 5.2, -73 


0.08 


Dec 2009 


1280 


32 


390 


2.8x2.1,53 


0.02 



Notes to Table[f]- a : Observed in dual 240/610 MHz mode. 
Table 2. Summary of the VLA archive observations 



Project Array Observation Frequency Bandwidth Integration FWHM, PA rms 

date (GHz) (IF1/IF2, MHz) time (min) (full array, "x", °) (/ijyb -1 ) 



AH406 


AnB 


Jul 1990 


8.41/8.46 


50 


40 


0.9x0.7,72 


20 


AT211 


A 


Apr 1998 


1.47/1.67 


25 


110 


1.3x1.1, -64 


25 


AT211 


B 


Aug 1998 


1.47/1.67 


25 


40 


4.6x4.1, -40 


40 


AT211 


A 


Apr 1998 


4.64/4.89 


50 


30 


0.4x0.4, -71 


50 


AT211 


B 


Aug 1998 


4.64/4.89 


50 


90 


1.4x1.4, -33 


25 


AT211 


C 


Apr 1998 


4.64/4.89 


50 


45 


7.2x4.1, 81 


40 



and instrumentation can be found in Jansen et al. (2001 and refer- 
ences therein). To provide a comparison with Chandra image anal- 
ysis, we reduced this data using SAS vl 1.0. 1 following the methods 
described in O'Sullivan et al. (2011). The EPIC instruments were 
operated in full frame mode, with the medium optical blocking fil- 
ter. Periods including background flaring, when the total count rate 
deviated from the mean by more than 3a, were excluded, leaving 
useful exposures of 12.0 ks (MOS1), 12.3 ks (MOS2) and 8.4 ks 
(pn). Point sources were identified using EDETECT_CHAIN, and re- 
gions corresponding to the 85 per cent encircled energy radius of 
each source (except that at the peak of the diffuse emission) were 
excluded. 

4 RADIO ANALYSIS 

The 1.28 GHz GMRT image at the resolution of 2.8" x 2.1" is 
presented in Fig. [2] as contours overlaid on the r'-band optical im- 
age from the Canada-France-Hawaii Telescope (CFHT). The im- 
age shows a double-lobe radio source with two fairly straight and 
weak jets and a dominant, compact component at the location of 
the cD galaxy. Fig.[3]presents the GMRT images at 610 MHz, 330 
MHz and 240 MHz. While the source is resolved into a triple at 
the 5"-resolution of the 610 MHz image, it is not spatially resolved 
in the images at 330 MHz and 240 MHz (FWHM~1 1" and ~14", 
respectively). The source has an angular size of ~ 27", correspond- 
ing to a linear size of ~150 kpc, fairly consistent at all frequencies 
and resolutions. 

Fig. [4^ shows the VLA image at 1.5 GHz, obtained from the 
combination of the A- and B-array data sets (Table[2j see also |Hines] 
|& Wills|1993> . The source morphology in this image is very sim- 
ilar to Fig. [2] The combined A+B-array contour image at 4.8 GHz 
is presented in Fig. [4]). Here, the brightest features are the central 
component and northern lobe; the S lobe is very faint and jets are 
barely visible. The central component and N lobe are also detected 
in the ~ l"-resolution image at 8.4 GHz (not shown here; Table|2]l. 

Fig.|5]zooms on the central region of CL09, by presenting the 
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Figure 2. GMRT contours at 1.28 GHz, overlaid on a CFHT r'-band im- 
age. The restoring beam is 2.8" X 2.1", p. a. 53°. The rms noise level is 
l<7=20^iJy. Contour levels are spaced by a factor of 2 starting from +3o\ 
Labels indicate the individual components of the radio source. 



VLBA image at 1.4 GHz. At VLB A resolution (8x6 mas), the cen- 
tral component is resolved into a quite symmetric double of ~35 
mas size (i.e., ~0.2 kpc). No clear core component is visible. The 
double accounts for a total flux density of 5.7 mJy, which is in 
agreement with 5.9±0.3 mJy measured for the central component 
in the VLA image at 1" -resolution (Table [5}. This implies that no 
further, extended structure is present in the central region in addi- 
tion to the VLBA double. 

The radio flux densities of CL09, obtained from the GMRT 
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contours start at +3cr and then scale by a factor of 2. Red dashed contours correspond to the — 3cr level. 



and VLA images presented above, are summarised in Table[3] along 
with the associated uncertainties. The flux density at 4.8 GHz was 
measured on the C-array image (not shown here; see Table [2](. At 
8.4 GHz we measured a lower limit on the flux, since only the cen- 
tral component and N lobe were detected. Inspection of the VLSS 
image at 74 MHz suggests the presence of a weak, unresolved 
source at the position of CL09. The source is not listed in the VLSS 
catalogue, implying that its flux density at 74 MHz is less then the 
5cj local rms level of the survey (1cj=100 mjy/beam, with FWHM= 
80" X 80"). We therefore assume an upper limit of 500 mjy for the 
flux density of CL09 at 74 MHz, also reported in Table [3] We also 
include two values from the literature; the 151 MHz measurement 
from t he 7C survey (|Riley et al.|1999fr , corrected to our adopted flux 
scale |Helmboldt et al. 2008 1, and a 15 GHz measurement made us- 
ing the Arcminute Microkelvin Imager (Hlavacek-Larrondo et al. 



2011). The 1.5 GHz radio power of the source, 1.1 x 10^ D W Hz" 1 
is consistent with the transition between FR I and FR II sources 
( |Bridle & Perley|1984[|Qw"en & Laing|1989l|Owen & White|1991| > 
as is often the case with optically luminous cluster and group- 
central radio galaxies. 



5 X-RAY ANALYSIS 
5.1 X-ray images 

We initially examined the structure of the ICM using the exposure 
corrected 0.3-3 keV image of CL09 shown in Fig. [6] As noted 
by |Iwasawa et al.| {200 1[ >, the ICM emission is brightest to the 
northeast, west and southwest of the core, with a depression to 
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Table 4. Radio properties of the individual components of CL09 



Component 


•Si. 28 GHz 


Sl.6 GHz 


Si. 8 GHz 


Ss^GHz 


«obs 


Ofit 




(mjy) 


(mJy) 


(mJy) 


(mJy) 






Central component 


7.06±0.35 


5.93±0.30 


1.51±0.08 


0.63±0.04 


1.28±0.04 


1.23±0.10 


Nlobe 


4.19±0.21 


3.57±0.18 


0.75±0.04 


0.38±0.03 


1.27±0.05 


, .,0 + 0.11 


S lobe 


1.55±0.08 


1.35±0.07 


0.46±0.03 




0.94±0.06 


90+ 20 


Njet 


3.45±0.22 


2.48±0.18 


0.46±0.03 




1.56±0.06 


1 44+0-06 


Sjet 


1.35±0.07 


1.17±0.06 


0.32±0.02 




1.12±0.06 


1 07+ - 20 



g 28.13 



o 
m 

Q 28.11 



i r 




□ 



09 13 45.496 45.495 



45.494 45.493 45.492 45.491 
RIGHT ASCENSION (J2000) 



Figure 5. VLBA image at 1.4 GHz of CL09. The restoring beam is 
0.0083" X 0.0063", p.a. -7°. Therms noise level is l(X=25/iJy. The peak 
flux density is 2.3 mJy beam -1 . Contour levels are 75 X [-1 (dashed), 1, 2, 
4, 8, 16] fiJy beam" 1 . 




Figure 6. Unbinned 0.3-3 keV exposure corrected Chandra image of the 
core of CL09, smoothed with a Gaussian of width 2 pixels (~1"). The 
black cross marks the optical centre of the BCG 



Table 3. Radio properties of NVSS J09 1345+405630 



Frequency, u 


S v (mJy) 


74 MHz 


<500 a 


151 MHz 


251±9 6 


240 MHz 


154.0±14 


327 MHz 


99.7±8.0 


610 MHz 


42.6±2.1 


1.28 GHz 


17.6±0.9 


1.5 GHz 


14.5±0.7 


4.8 GHz 


3.5±0.2 C 


8.4 GHz 


>l d 


15 GHz 


0.80±0.04 e 


Q 151 MHz-15 GHz 


1.25±0.01 



Notes t o Table [3]- a : 5<t upper limit from the VLSS; b : from the 7C 
survey ( Riley et al. 1999), converted to our adopted flux density scale using 
the conversion factor listed by Helmboldt et al. 1 2008 1; c : measured on the 
C-array image; d : lower limit since only the core and N lobe are detected 
in the 8.4 GHz image; e : Arcminute Microkelvin Imager measurement 
reported by Hlavacek-Larrond o et al.|j2011} . 



the northwest. This hole was first observed in ROSAT HRI data 
( |Fabian & Crawford|1995l|Crawford & Vanderriest|l"996l l. The hole 
is bounded on its northwest side by a linear filamentary structure. 
To determine the significance of the filament, we placed ~ 1 .5" x 7" 
box regions on the filament and immediately to either side, and de- 
termined the number of counts in each region. The exposure cor- 
rected 0.3-3 keV surface brightness of the filament is 2a greater 
than the region inside (in the hole) and 3.6<r greater than the sur- 
face brightness immediately outside. 

As well as the northwest hole, HL11 identify a second cav- 
ity to the southeast. To provide a clearer view of the structures 
in the ICM, we modelled the surface brightness distribution of 
the cluster in the CI AO SHERPA package using an elliptical /?- 
model, with an additional constant component to represent the 
background. A monoenergetic exposure map with energy 1.39 keV 
(chosen to match the mean photon energy of the 0.3-3 keV im- 
age) was used to correct for exposure variations across the source. 
The effect of the PSF is probably minimal and was therefore 
neglected. The best fitting model had core radius 4.42"±0. 11 
(24.3±0.6 kpc), /3=0.578±0.004, ellipticity 0.136±0.010, and po- 
sition angle 9.3°±2.2, indicating that the major axis of the ICM 
runs almost east-west. All parameters of the /3-model were allowed 
to vary freely in the fit, including the central position. Inclusion of a 
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Figure 7. Chandra 0.3-3 keV residual image after subtraction of the best fit- 
ting elliptical /3-model, smoothed with a Gaussian of width 2 pixels (~1"). 
The northwest filament and optical galaxy centre are labelled as in Fig. [6] 
Dashed cyan circles or ellipses indicate regions used to estimate cavity 
properties. 1.28 GHz GMRT contours are overlaid (same as Fig.[2). 



second /3-model had no significant effect on the model parameters, 
the second component primarily modelling emission in the core of 
the BCG. 

Fig. ^] shows the residual image created by subtracting the 
best-fitting model from the 0.3-3 keV image, smoothed with a 2"- 
width Gaussian to bring out small-scale features. A central bright 
source is visible coincident with the BCG core and inner compo- 
nent of the radio source. There is a bright extension out to ~5.5" 
northeast, almost perpendicular to the radio jet axis, which is corre- 
lated with the brightest of the gaseous filaments north of the BCG 
and with a small spur of 1.28 GHz radio emission (see Fig. Wt. 
However the strongest features are two dark regions of negative 
residuals to northwest and southeast, roughly coincident with the 
jet axis, but not extending as far as the radio lobes. These are the 
cavities, where thermal ICM plasma has presumably been driven 
out by the relativistic plasma of the radio jets, leaving a deficit in 
X-ray emission. The northwestern cavity corresponds to the hole 
seen in the raw image, and has bright residuals extending around its 
western and northern sides, including the northwest filament. These 
bright features may arise from ICM gas pushed out of the cavity. 
The southeast cavity is less well-correlated with the jet structure, 
extending to the east beyond the radio emission. There is also a 
weak negative residual to the southwest of the central component, 
which is uncorrelated with the jet. Our estimates of cavity size and 
position differ somewhat from those of HL1 1. We discuss these dif- 
ferences (and the physical properties of the cavities) in Section [73] 

Heavier smoothing reveals residual features on larger scales. 
Fig.|9]shows the same Chandra 0.3-3 keV image smoothed with a 
~6"Gaussian. The northwest cavity, surrounded by brighter emis- 
sion, is still visible as a dark hole in the image centre. However, 
the heavier smoothing brings out negative residuals running from 
the north around the west side of the core and connecting to the 
west side of the southwest cavity. The weak negative residuals seen 
in Fig. [7] southwest of the core and extending east of the southeast 




Figure 8. HST 622W image of the core of CL09, overlaid with 1.28 GHz 
radio contours (solid red, as in Fig.[2j and contours of smoothed positive 
X-ray residuals (dashed black, starting at 0.5 counts/pixel and increasing in 
steps of factor 2). 

cavity are thus likely to be part of a larger structure, rather than cav- 
ities. The positive residual to the northeast of the core is also found 
to extend further, curving around the eastern and southern side of 
the core. 

To check the robustness of the features revealed by the surface 
brightness modelling, we performed a separate fit to a 0.5-3 keV 
XMM-Newton image (Fig. [9}. The northwest cavity is not resolved, 
as its ~4" radius is smaller than the XMM point-spread function. 
However, the spiral structures seen in the Chandra image are clear, 
and can be traced in both positive and negative residuals curving 
out from the core to ~300 kpc. 

Measuring the significance of the cavities and spiral structure 
is difficult, since they affect the ICM as far as it can be clearly 
traced in the images. We therefore compare radial profiles across 
the features with the best-fitting surface brightness model. Fig. 1 10| 
shows profiles in four directions, two along the radio jets, and two 
roughly perpendicular to them. Surface brightness deficits are visi- 
ble at radii of 2-6" along the northwest jet (the hole) and 5-7" along 
the southeast jet. In individual bins, the deficits are significant at the 
4-5cr level. The northwest filament is visible as a surface brightness 
excess at 7-10" radius. Neither profile shows a significant deficit at 
the position of the radio lobes. A significant excess is seen at radii 
<5" in the northeast profile, corresponding to the residual emis- 
sion extending in this direction from the core. The western profile 
has an excess at 5-10", agrees with the average surface brightness 
profile at 10-20", and rises above it at larger radii. This corresponds 
to the residual features observed; first a surface brightness excess 
on the west side of the northwest cavity, then the region of low rel- 
ative surface brightness which forms part of the negative residual 
spiral in Fig. [9] These results confirm that the cavities and associ- 
ated small-scale structures are both real and significant, and provide 
some support for the spiral features observed in the residual maps. 

5.2 Gas properties 

To determine the physical properties of the ICM, we extracted spec- 
tra in radial bins with widths chosen to ensure at least 3000 net 
source counts per annulus. A central 1" radius region was excluded 
to remove emission from the AGN. The spectra were fitted with an 
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Figure 9. Upper panels: Chandra 0.3-3 keV {upper left) and XMM-Newton 0.5-3 keV (upper right) residual images after subtraction of the best fitting elliptical 
/3-model and Gaussian smoothing. The Chandra smoothing scale is 12 pixels (~6")> the XMM-Newton scale in 3 pixels (13.2"). The two images have the 
same angular scale and alignment, and the green cross in the Chandra image indicates the optical centre of the BCG. 1.28 GHz GMRT contours are overlaid 
in green on the centre of the XMM-Newton image. The same spiral pattern of residuals is visible in both images. The surface brightness hole to the northwest 
of the core is also visible in the Chandra image, but is too small to be resolved in the XMM image. Lower panels: 0.3-3 keV Chandra residual maps, smoothed 
with Gaussians of width 1" (lower left) and 6" (lower right). Cyan contours indicate the 1.28 GHz radio emission. Green partial annuli indicate regions used 
to create the radial profiles shown in Fig. |10| 



absorbed APEC thermal plasma model ( |Smith et al.|200l) , using 
the XSPEC projct model to account for the effects of projection. 
Abundances were tied between some bins to reduce uncertainties 
and the typical value was found to be ^0.6Zq, declining outside 
~40" (~220 kpc). Fig. 11 shows radial profiles of the various gas 
parameters derived from the fit. Our temperature, electron density 
and entropy profiles are in good agreement with those of HL11. 
We are able to trace the ICM to ~150" (~850 kpc) within which 
the cluster has a bolometric X-ray luminosity ~2xl0 45 erg s _1 . 
The radio jets extend to ~14" (77 kpc), corresponding roughly to 
the outer radius of the third annulus. The luminosity within this ra- 
dius is -~8.5x 10 44 erg s _1 . The peak temperature is 7.7+g'g keV, 



falling to 3.9±0.3 keV in the central 5" (~28 kpc), where the cool- 
ing time is ~ 1 Gyr. 

To check for any indication of strong heating or cooling in the 
central region, we also fit projected APEC models to spectra ex- 
tracted from four l"-wide radial bins, starting at 1.5". We find a 
close agreement between the projected and deprojected core tem- 
peratures, with no indication of a change in temperature gradient 
from that observed in the deprojected profile. In general the tem- 
perature profile shows that the heating the AGN has delivered (ei- 
ther radiatively or mechanically via the jets) has been insufficient 
to disrupt or destroy the strong cool core of the cluster. 

We cannot accurately estimate the rate at which gas is cool- 
ing out of the ICM, since the XMM RGS data are not deep enough 
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to allow high-resolution spectral modelling. We can estimate mass 
deposition rates based on the Chandra ACIS spectra, but these are 
very likely to be overestimated (since the calculation neglects heat- 
ing and other important physical effects) and should be considered 
upper limits. Replacing the single temperature APEC model for the 
central bin in the deprojection analysis with a cooling flow model 
(MKCFLOW, |Mushotzky & Szymkowiak|1988J we find a mass de- 
position rate M=319+34 M§ yr _1 within ~27 kpc. However, the 
fit is not a significant improvement over the APEC model and has 
an unrealistic maximum temperature (~11 keV). Alternatively, we 
can estimate the isobaric cooling rate: 



M 



2fj,m p L) 



5kT 



(1) 



where y, is the mean molecular weight (0.593), m p is the pro- 
ton mass, Ljf.boi is the bolometric X-ray luminosity of the central 
spectral bin, and kT the corresponding temperature. This gives an 
estimate of M ~235 M Q yr" 1 . 

Using the mapping technique described in O'Sul livan et al.| 
< [20 1 1 1 >, we created spectral maps of the cluster core. Each 
2.5x2.5" pixel of the map represents an absorbed APEC model 
fitted to a spectrum extracted from a circular region chosen to con- 
tain at least 1000 net counts. The extraction regions had radii of 
2-30" and therefore the spectra for neighbouring regions are gen- 
erally not independent. Fig.|12|shows the temperature map. Uncer- 
tainties on fitted temperatures are 7-17%. The map shows that the 
coolest part of the ICM is extended along a northeast-southwest 
axis, roughly perpendicular to the (projected) axis of the radio jets, 
and not aligned with the major axis of the best fitting surface bright- 
ness /3-model. However, we note a good correlation between cool 
temperatures and optical filaments; all the filaments are found in 
regions with kT<5 keV, and the complex of filaments northeast of 
the BCG extends along the same axis as the coolest part of the ICM. 

The spectral map does not show any structures correlated with 



Table 5. Characteristic total masses, radii and gas fractions for the cluster. 
Overdensity radii are derived from the X-ray mass model, X-ray (X) and 
weak lensing (wl) masses are measured within those radii and gas fractions 
are derived using the X-ray mass estimates. 



Overdensity, A 



Ra 

(kpc) 



Ma.x 
10 14 M, 







10 14 M 



•-gas 



O 



2500 

500 

200 



487±29 2.60±0.50 
1088±97 5.83±1.69 
1648±164 8.10±2.64 



2.14±0.43 0.09±0.01 
4.12±0.98 0.11 ±0.02 
0.12±0.03 



the cavities, northwest filament, or spiral structures discussed in 
section [5~T| Extracting spectra from the cavity regions shown in 
Fig. [7] and from regions chosen to cover the northwest filament 
and residuals around the northern cavity, we found no signifi- 
cant temperature differences beyond those expected from the map. 
Similarly, placing regions to match the spiral residuals shown in 
Fig. [5] we find only a 1.5(7 significant temperature difference be- 
tween the positive and negative residuals, the negative region being 
marginally cooler (5.78±0.30 keV compared to 6.35±°' J4 keV). 



5.3 Mass Profiles 

We determine the total gravitational mass profile of the cluster, and 
its gas mass profile, using the Joint Analysis of Cluster Observa- 
tions system (JACO, |Mahdavi et al.||2007) . This allows us to fit 
models jointly to the Chandra and XMM-Newton data, taking ad- 
vantage of the improved signal-to-noise ratio thus available. The 
dark matter distribution was modelled as an NFW (Navarro et al. 
1 1995 1 [Navarro et al.|1997] > profile with a freely fitted concentration 
parameter. For the gas profile, we used a model consisting of three 
/3-models, with the central component multiplied by a power law 
with free slope to model the bright emission from the central cool 
core. The minimum radius to which the mass profile can be deter- 
mined is limited by the point-spread function (PSF) of the XMM- 
Newton EPIC instruments. Since the angular size of the BCG is 
small (~3" across), we cannot resolve the mass profile within the 
galaxy. We therefore do not include an additional mass model to 
account for the stars, but consider them to be included in the NFW 
model. 

The results of this fitting are shown in Fig. [T3] For compari- 
son, we have also included mass data points derived from the tem- 
perature and density profiles shown in Fig.[TT] These are calculated 
under the assumption of hydrostatic equilibrium, as is the case with 
the JACO modelling, but using the apparent gradients between bins. 
The two estimates agree reasonably well. From the fitted model we 
are able to estimate characteristic masses and radii for CL09, as 
well as the enclosed gas fraction, for different degrees of overden- 
sity. These are shown, with lcr uncertainties, in table [5] The gas 
fraction within R2500 is well within the range of gas fractions found 
for the sample of massive clusters studied by |Allen et al.| ([2008 ). 

We also include a comparison with weak lensing mass es- 
timates drawn from the Canadian Cluster Comparison Project 
(CCCP). This is a multiwavelength study of a sample of ~ 50 X- 
ray luminous clusters of galaxies, which includes CL09. To deter- 
mine weak lensing masses, deep imaging data were obtained us- 
ing Megacam on the Canada-France-Hawaii Telescope. The weak 
lensing analysis is described in Hoekstra et al. (in prep). A detailed 
discussion of the analysis and the various issues that arise when in- 
terpreting the data is provided in Hoekstra (2007 1, which analysed 
an initial sample of 20 clusters. 
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Figure 11. Radial profiles of gas temperature, electron number density (n e ), metal abundance, pressure (P), entropy (S) and isobaric cooling time (T coo j) for 
the cluster. Black points represent the deprojected fits. For comparison, projected temperatures in four bins at 1 .5-5.5" are marked in grey. Errorbars or regions 
indicate \u uncertainties. 



Fitting the measurements between 0.5-2 Mpc with a singu- 
lar isothermal sphere model, we find a best fit which corresponds 
to a velocity dispersion of <j=974^}gQ km s~ 1 . We also compute 
aperture masses, which we convert to spherical masses under the 
assumption that the density profile can be described by an NFW 



profile (see Hoekstra 2007 for details). The resulting masses are 
marked on Fig. 1 13| Table[5llists characteristic masses measured at 
values of ta determined from the X-ray analysis (note we cannot 
determine A/200 this way). The uncertainties include the statistical 
error caused by the intrinsic shapes of the sources and the noise 
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Figure 12. Temperature map of CL09, in units of keV. Uncertainties are 
7-17%, with the central regions having smallest uncertainties. 1.28 GHz 
GMRT contours are overlaid. 
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Figure 13. JACO profiles of total gravitational mass (solid line) and gas 
mass (dashed line) for CL09, extending to R500 ■ Grey regions indicate 1 a 
uncertainties, though for gas mass these are similar to the width of the line. 
The profiles are derived from a joint fit to the Chandra and XMM-Newton 
data, under the assumption of hydrostatic equilibrium. For comparison, we 
include mass points derived directly from our temperature and density pro- 
files (diamonds) and from a weak lensing analysis (errorbars, Hoekstra et 
al. in prep.). 
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Figure 14. Integrated radio spectrum of CL09 between 74 MHz and 4.8 
GHz. The solid line is the CI fit to the data. Open points and limits indicate 
points included in the fit. Two points deviant points not included in the fit 
are marked by solid symbols. 

from distant large scale structure (e.g. Hoekstra et al. 2011). The 
lensing results are in good agreement with the X-ray mass determi- 
nations. 



6 RADIO SPECTRAL ANALYSIS AND PHYSICAL 
PARAMETERS 

6.1 The integrated spectrum 

In Fig. [14] we show the integrated radio spectrum of CL09 be- 
tween 74 MHz and 15 GHz, obtained using the flux densities 
in Table [5] The source has a straight spectrum with steep slope, 
a = 1.25 ± 0.01, between 151 MHz and 15 GHz. The upper 
limit at 74 MHz suggests a low-frequency flattening of the spec- 
trum, with a ^ 1 below 151 MHz. We fitted the spectrum using a 
continuous injection (CI) model (e.g., Kardashev 19621, in which 
electrons are continuously injected in a region of constant magnetic 
field. The CI best-fit model (solid line in Fig. provides an in- 
jection spectral index aiuj = 0.75 ± 0.05 and a break frequency 
^brcak — 51^3° 6 MHz. Two additional data points are marked on 
Fig. [141 at 327 MHz (from the WENSS survey, |Rengelink et al.| 
|1997| > and 28.5 GHz ( |Coble et al.|[2007] >. It is unclear why these 
measurements disagree with the other data points, and we there- 
fore do not include them in the fit. 

We note that a spectral break at a frequency higher than 15 
GHz is implausible for energetic reasons. Assuming 5=2ainj+l, 
where 5 is the slope of the injection energy spectrum of the elec- 
trons, and vtreak ^ 15 GHz, the observed spectral index at lower 
frequency (a = 1.25 « oimj) gives 5=3.5, which implies an ex- 
tremely steep spectrum of the injected electrons and an untenably 
large energy budget of the electron population, contributed by elec- 
trons at lower energies. 

Fig. [T5] shows the integrated spectra of the individual compo- 
nents of CL09, computed from the flux densities in Table|4] Given 
the absence of an obvious spectral curvature and limited number of 
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data points, we used a simple power-law model to fit the spectrum 
of each component. The slopes provided by the best-fit models are 
summarised in Table[4] The spectral indices of the individual com- 
ponents are all steep (a > 1), including the central component. 



6.2 Physical parameters 

We calculated the physical parameters of the source using the stan- 
dard assumptions that the relativistic particle and magnetic field en- 
ergy densities are uniformly distributed over the radio source vol- 
ume and in approximate energy equipartition. We also imposed a 
low-energy cutoff 7mm (where 7 is the electron Lorentz factor) in 
the energy distribution of the radiative electrons, instead of adopt- 
ing a fixed interval of frequency (typically 10 MHz-100 GHz; e.g., 
|Pachol czyk 1970). It has been shown that a frequency cutoff of 10 
MHz may neglect the input of 7 < 10 3 — 10 4 electrons for typical 
equipartition magnetic fields of ~ 1—10 /iG (Brunetti et al. 1997 
Beck & Krause 2005 1. Since a large fraction of the relativistic par- 
ticle energy may be at lower 7, imposing a low-energy cutoff allows 
to take into account the contribution from these particles. In our cal- 
culations, we adopt 7mm=100, which corresponds to a low-energy 
cutoff of ~50 MeV. 

We adopted an injection spectral index a m j=0.75 provided by 
the CI best-fit model (see Section [6T| l, and normalised the model 
to the radio luminosity at 151 MHz, the lowest frequency at which 
the source is detected, and thus the least affected by spectral age- 
ing of the electron population. Cylindrical geometry was adopted 
to calculate the volume of the source. The results are summarised 
in Table [6] where we list the radio luminosity at 151 MHz, injec- 
tion spectral index ai n j , break frequency, volume, and equipartition 
magnetic field B aq . The last three columns provide estimates of the 
radiative age f ra d, velocity growth « gr owth/c and integrated lumi- 
nosity of the source (see Sect.|6.3}, 



6.3 Radiative age 

Under a number of assumptions, the knowledge of the break fre- 
quency in the spectrum of a radio source allows us to estimate the 
time elapsed since the source formation (e.g., Myers & Spangler 
|1985[ >. If radiative losses dominate over expansion losses, the mag- 
netic field is uniform across the source and remain constant over the 
source life time, and reacceleration processes can be neglected, the 
radiative age (the time since the AGN last powered the jet, launch- 
ing relativistic electrons into it) can be obtained as 



= 1590 



(B 2 + B 2 



-[(1 + z)i/ b reak] 



Myr 



CMB / 



where fbrcak is expressed in GHz, and B cq and Bcmb in /j,G 
(Parma et al. 2007}. -B eq is the equipartition magnetic field, and 
Bcmb = 3.2(1 + z) 2 is the equivalent magnetic field of the cosmic 
microwave background (CMB) radiation, i.e., the magnetic field 
strength with energy density equal to that of the CMB at the red- 
shift z. 

Using the break frequency in Table 5, we estimate a radiative 
lifetime of ~ 130 Myr for CL09. We also derive an estimate of the 
source growth velocity v gr owth/c, where c is the speed of light, ob- 
tained assuming a constant velocity and a linear size of LLS=150 
kpc, as measured from the radio images in Sect. [4] i.e., u gr0 wth = 
LLS/trad- We found u growth < 0.004c (~1200 km s" 1 ). For 
comparison, the ICM sound speed varies from ~870-l 120 km s~ 
over the length of the jets. 



7 SUMMARY AND ANALYSIS OF RESULTS 
7.1 The cluster 

Our X-ray observations show CL09 to be a luminous, massive clus- 
ter with a strong cool core. The ICM is only mildly elliptical in the 
plane of the sky, with no visible subclumps or near neighbours. 
Spectral mapping shows a degree of ellipticity in the temperature 
distribution, particularly in the core. To some extent this ellipticity 
may be a product of the jet/ICM interactions, with the inflation of 
cavities reducing the amount of cool gas along the jet axis. How- 
ever, there is no indication of azimuthal asymmetries outside the 
core which might invalidate our temperature profile or mass anal- 
ysis. In general, it seems fair to describe the cluster as essentially 
relaxed, despite the structures detected in the ICM. 

In order to look for more subtle deviations, we compared the 
physical properties of CL09 to other clusters via the M-T and L-T 
relations. We took the L-T relation for the REXCESS cluster sur- 
vey ( |Pratt et al.|20 09 1, specifically the relation between temperature 
and luminosity in the 0.15-1 XR500 radial range. Using the value of 
R500 determined in Section [53| we extracted a Chandra spectrum 
and fitted an absorbed APEC model. The resulting luminosity falls 
almost exactly on the relation. Following the same approach, we 
compare with the M-T relations of Vikhlinin et al. (2006} and |Ar-| 
|naud et al.| l |2005) . In both cases the cluster falls on the relation to 
within the la uncertainties. 

The weak spiral structure found in the residual images may be 
an indicator of gas motions in the ICM. Such structures are typi- 
cal of sloshing motions, where the cluster core is disturbed by the 
nearby passage of another massive object ( Ascasibar & Mar kevitch] 
|2006| |ZuHone et al.|20lT) . Spiral structures are observed in both 
surface brightness and temperature in a number of disturbed galaxy 
clusters and groups (e.g., |Clarke et al.|2004{ |Randall et al.|2009| 



|Lagana et al . 2010, Blanto net al.|20lT i. We do not observe the 
expected spiral structure in the temperature map of CL09, nor do 
we observe the sharp paired cold fronts observed in many sloshing 
systems. However, this may simply be because the Chandra ob- 
servation has insufficient depth to allow fine-grained temperature 
mapping. 

There is no clear evidence of a perturber in the X-ray data, 
but we can consider what kind of system might be required to dis- 
turb the cluster. Disturbance will be more easily caused by higher 
mass systems, but given the mass of CL09, M2oo=8-1 x 10 14 M Q , 
it seems likely that a group-scale system with a mass of a few 
10 13 Mq would be required. Such a group would be expected to 
have a luminosity > 10 43 erg s _1 , and would thus be detected if it 
were still close to CL09. As no perturber is detected, we are left 
with the possibility that either the disturbing system had a very low 
mass (e.g., an individual infalling galaxy or poor group), that the 
gas associated with the disturber has been stripped, that it is out- 
side the ACIS-I field of view, or that it has already fully merged 
with CL09 leaving only the spiral residuals as a record. Simula- 
tions of cluster-cluster mergers also show spiral structures in cases 
where the merger is off-axis (e.g., Poole et al. 2006). The cool core 
and weak spiral pattern in CL09 suggest that if a merger took place, 
the cluster cores coalesced some time ago and the cluster is now 
largely relaxed. 



7.2 The radio source 

Our analysis of the GMRT and VLA observations of the large-scale 
radio source provide a relatively clear picture of its status. The 
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Figure 15. Integrated radio spectra of the individual components of CL09. The solid line in each panel is the power-law fit to the data. The image is the VIA 
1.5 GHz image (same as Fig.[4j(). 



Table 6. Physical parameters of CL09 
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source is an FR I/FR II transition system, with relatively straight, 
resolved jets and fairly small, bright lobes. The source has a steep 
spectral index between 151 MHz and 15 GHz (q=1.25±0.01), 
with a probable break at ~50 MHz. Separate spectral indices mea- 
sured in the jets, lobes and central component are also steep and 
broadly comparable, with no clear change of spectral index with 
distance from the centre. The source is therefore probably old 



(t rad =130 Myr) and no longer powered by the AGN. There is no 
indication of faint, low-frequency emission on scales larger than 
the jets and lobes, and therefore no indication of any earlier radio 
activity. 

The VLBA observation reveals a small-scale (~200 pc) double 
source, both components of which are extended. No point-like core 
is detected, indicating that both components of the double source 
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are probably lobes, rather than being associated directly with the 
AGN. The flux from the VLBA double is comparable to the flux 
of the central component of the large-scale source, with no indi- 
cation of any additional emission on other scales. The similarity in 
flux suggests, firstly, that the inner double is a separate structure 
from the large-scale jets. There is no indication of intermediate- 
scale structures between the inner double and the base of the large- 
scale jets 20-30 kpc from the centre. The well-defined, compact 
lobes of the inner double also support this conclusion. Secondly, 
if the central component of the large-scale source is merely the 
unresolved inner double, the steep spectral index measured for the 
central component is also applicable to the VLBA double, implying 
that it too is old. Double sources with such small physical scales 
are often assumed to be young, but it has been suggested that many 
such sources may be short-lived ( Kunert-Bajraszewska et al.|2010| >, 
and examples where spectral index measurements show that the jets 
have already shut down are not unknown (e.g., |Orienti et al.|2010| >. 
The age of the inner double cannot be confirmed without new VLBA 
observations at other frequencies. Such observations would also de- 
termine whether the central AGN is truly radio quiet or whether it is 
self-absorbed at 1.4 GHz, which is not uncommon for radio cores. 

In the plane of the sky, the axis of the VLBA inner double 
differs from that of the large-scale jets by ~17°. The similar sizes 
of the two components of the inner double suggests that it is not 
strongly relativistically beamed. Under the assumption that the jets 
are intrinsically symmetric, the jet-to-counterjet brightness ratio 
can be used to constrain the velocity /3c (where c is the speed of 
light) and orientation 8 of the jets (see Giovanni ni et al.|1994| and 
references therein). We measure a ratio R = 1.3 between the N 
and S lobe flux densities of the double. Adopting a = 1.23, i.e., 
the spectral index of the unresolved central component on arcsec- 
ond scale (Tab.4), we estimate f3cos8=0.03. Given the distribution 
of jet velocities on the milliarcsecond scale, where usually B > 0.5 
(e.g., |Giovannini et al.|200i] >, it is very likely that the inner double 
is lying almost in the plane of the sky, as 9 is larger than 80° for 
8 > 0.2. 

7.3 Cavity energetics and timescales 

From the surface brightness analysis, we have fairly strong evi- 
dence of a cavity corresponding to the northern jet and another 
whose position coincides with the southern jet but which is some- 
what larger than the current radio emission. The southern cavity 
lies within the negative spiral residual, and it is therefore difficult 
to reliably determine its size. Neither of these cavities is correlated 
with the radio lobes, only with the jets. 

We might also have expected evidence of cavities coincident 
with the radio lobes. Although no clear X-ray holes are seen at the 
positions of the lobes, we can test whether weak surface bright- 
ness depressions might be present. Adopting the regions shown in 
Fig. [7] (circles of radius 13.8 kpc) we compare the exposure cor- 
rected surface brightness with equally sized regions placed to ei- 
ther side of the lobes, along an ellipse defined by the best fitting 
surface brightness model. The southern lobe region has a surface 
brightness within la of its neighbours. The northern lobe region 
is actually brighter than the neighbouring regions, since it contains 
part of the northeast filament. We therefore have no evidence of 
cavities coincident with the radio lobes. 

We can also test whether our data are sufficiently sensitive to 
detect cavities of the expected size at the position of the lobes. We 
estimate the flux deficit expected for cavities of size comparable 
to the radio lobes, based on the mean density and temperature at 



that radius. We predict that such cavities should produce a deficit 
of ~30 counts (0.3-3 keV), which is equivalent to only a 2.9a drop 
for the southern lobe. We cannot therefore be certain that we would 
reliably detect cavities, if they do exist. A comparable calculation 
for the northern cavity (i.e., the X-ray hole) suggests that we should 
expect a decrement of ~400 counts. There are no undisturbed re- 
gions at that radius which can be used for comparison, but the ob- 
served flux is ~320 counts lower than that predicted by the best- 
fitting surface brightness model. 

Inverse-Compton scattering of cosmic microwave back- 
ground (CMB) photons by electrons in the relativistic plasma of 
the radio lobes could potentially produce additional X-ray flux in 
the areas where we expect to see cavities. Under the assumption 
of equipartition, we estimate the energy of the electron population 
of each lobe to be 5-10x 10 56 erg, for an injection spectral index 
a irlJ =0.7, minimum Lorentz factor 7 m in=100, and an equal divi- 
sion of the energy in the relativistic particle population between 
electrons and protons or positrons (k=l). Following the method 
described in O' Sulli van et al.| p010), we estimate the expected 
inverse-Compton flux from this electron population using the re- 
lations of Erlund et al. (2006). The expected 0.3-3 keV flux from 
each lobe is 1.3-2.2X 10 _1Y erg s _1 cm -2 , a factor ~10~ 4 lower 
than the flux measured in these regions. It is therefore clear that our 
failure to detect cavities is not caused by inverse-Compton emis- 
sion. We note that the synchrotron self-Compton emission from 
the radio lobes will be even fainter, since the energy density of ra- 
dio photons in the lobes is a factor ~10~ 3 lower than that of the 
CMB. 

Since we cannot be certain that the radio lobes do not have 
associated cavities, we estimated the energy required to inflate the 
observed cavities, and any potential cavities at the position of the 
lobes. These estimates are shown in Table|7] as are dynamical esti- 
mates of the timescale over which the cavities have formed and the 
implied mechanical power of the radio jets. 

Cavity volumes and the uncertainty on the volume are esti- 
mated following the methods described in O'Sul livan et al.| ( [201 1| >. 
We assume the jets to be in the plane of the sky, and parameters 
from our radial spectral fits, which account for projection effects. 
We define the energy of each cavity to be 4pV, the enthalpy of a 
cavity filled with a relativistic plasma. We estimate the time for the 
cavity to rise to its current position at the speed of sound (t son ic). 
the time for it to buoyantly rise to its position (tbuoy, |Churazov| 
|et al.||200T) , and the time for the ICM to refill the current cav- 
ity volume (t re /izz). For comparison, the bolometric X-ray lumi- 
nosity within the radius to which the radio jets extend (77 kpc) is 
~8.5x 10 44 erg s _1 . The mechanical energy required to inflate the 
observed cavities (pV) is greater than this luminosity, assuming the 
commonly used buoyant timescale. Even if we use the longer radia- 
tive age estimate (~ 130 Myr), the enthalpy of the cavities (4pV) 
still exceeds the cooling luminosity. The power of the jets is there- 
fore more than sufficient to balance cooling in the cluster core. If 
we assume that all the energy of the cavities is available to heat the 
ICM, the jets can balance cooling out to at least 200 kpc (for the 
refill timescale), or out to the edge of our temperature profile (for 
the buoyant timescale). 

Our estimates of the cavity sizes and positions differ from 
those of HL11. We find the northwest cavity to be somewhat 
smaller (HL1 1 use a cavity region which extends across the north- 
west X-ray filament, which we consider likely to be the cavity rim) 
and the southeast cavity to be of similar size but closer to the cluster 
core. These differences arise from our different approaches to iden- 
tifying cavities. HL11 examine unsharp-masked and residual im- 
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Table 7. Cavity sizes, powers and timescales 

Cavity radius D core 4pV t sonic t buoy t re 

fill Psonic Pfmoy ^refill 
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Figure 16. Cavity power vs. integrated 10 MHz-10 GHz radio power. Sys- 
tems in the sample of |Birzan et al.| j2008) are marked by grey triangles, 
those in the sample of |O ; Su7h'v1metal.|j201 1) by black circles. CL09 is 
marked by blue squares, the filled square indicating the radio power ob- 
served, the open square the radio power after correcting for radiative age- 
ing, lcr uncertainties on radio power and cavity power (calculated using 
the buoyancy timescale) are indicated by the error bars. For the O'Sullivan 
et al. groups and CL09 additional narrow-width error bars indicate the la- 
error range allowing for alternate measures of cavity age (sonic and refill 
timescales). The dotted and solid lines show the best-fitting relationships 
found by |Brrzan et"aT](2008) and |Q' Sullivan et al.|(20"TT) respectively. 



ages, but select the cavity regions from the raw image, whereas we 
select regions from the residual image. This demonstrates the inher- 
ent subjectivity of cavity identification. HL11 estimate timescales 
and powers for the cavities which also differ from ours, but our es- 
timates of total jet power are in reasonable agreement, given the un- 
certainties. One potentially important difference in interpretation is 
that the outer edges of the cavity regions used in HL1 1 overlap the 
radio hotspots, and they therefore say that the cavities and hotspots 
coincide. Our regions do not overlap the radio lobes, and both their 
cavity regions and ours are poorly correlated with the radio lobes. 

HL1 1 also compare CL09 to the known relations between jet 
power, radio luminosity and cooling luminosity. We perform the 
same comparison, using our estimates of these parameters. Fig.|16| 
shows a plot of cavity power (for the detected cavities only) against 
integrated 10 MHz-10 GHz radio luminosity, which we estimate 
to be 2.1 x 10 42 erg s -1 . CL09 falls above the best fitting relation, 
but is comparable to a number of other clusters with large jet pow- 
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Figure 17. Cavity power (4PV) against cooling luminosity (bolomet- 
ric X-ray luminosity of the volume within which the cooling time 
tcool ^7.7 Gyr). The blue square indicates CL09, using the power esti- 
mated from the two detected cavities. Black crosses indicate the elliptical 
galaxies and groups sample of |Cavagnolo et aL] (2010} and grey triangles 
the cluster sample of Birzan et al. (2008). lcr uncertainties on cavity power, 
calculated using the buoyancy timescale, are indicated by vertical errorbars. 
lcr uncertainties on L coo ; are indicated by horizontal bars for CL09 and 
the Birzan et al. sample, but in many cases the statistical uncertainties are 
smaller than the data points. Cavagnolo et al. do not quote uncertainties on 
cooling luminosity. For CL09, additional narrow-width bars indicate the lcr 
uncertainties on cavity power allowing for alternate measures of cavity age 
(refill or sonic timescales) and the change in L coo ; between the neighbour- 
ing bins of the spectral profile. The solid line indicates equality between 
L coo j and V cav , while the dotted lines indicate equality if PV or 16PV of 
energy were required to balance cooling. 

ers. Its offset from the relation is in part due to spectral ageing 
which has both steepened the spectral index and lowered the flux 
we measure across all frequencies. We can partially correct for this 
effect by estimating the radio luminosity using the injection spec- 
tral index ai n j=0.15 normalised to the 151 MHz flux, which (of the 
frequencies at which we detect the source) is least affected by radia- 
tive ageing. We estimate this partially corrected radio luminosity to 
be 3.6x 10 42 erg s~ . This moves the source closer to the best fit- 
ting relation, but only by a small amount. A more important issue 
is the timescale used. All data points in Fig. [T6] use the buoyancy 
timescale, but for CL09 a longer timescale may be more appropri- 
ate, given the longer radiative age estimated from the radio spectral 
analysis, and the possible link between the onset of star-formation 



© 2012 RAS, MNRAS OOO.fTp?] 



A type 2 QSO in a cooling flow 17 



and shut-down of the jets (discussed in Section |8T) . This would 
move the cluster downward, toward the best-fitting relation. How- 
ever, we note that if there are undetected cavities coincident with 
the radio lobes, including them would move the point further away 
from the relation. 

Fig. [IT] shows a plot of cavity power against L COD i, the lumi- 
nosity of the region within which the cooling time is <7.7 Gyr. 
CL09 has a cavity enthalpy in excess of the cooling luminos- 
ity, falling on the PV line, indicating that the mechanical power 
needed to inflate the cavities is sufficient to balance cooling with- 
out any contribution from the energy stored in the relativistic par- 
ticles or magnetic field of the radio jets. CL09 is quite comparable 
to the clusters of the|Birzan et al. (2008 1 sample, falling within the 
scatter of the cluster population. In both Fig.[l6]and[T7j our results 
are qualitatively similar to the findings of HL1 1, suggesting that the 
radio source in CL09 is capable of balancing radiative cooling and 
that it is mildly underluminous for its jet power. 

It should be noted that the jets are very likely to have expanded 
supersonically for some portion of their history and would have 
driven shocks into the ICM during this period. Shocks associated 
with radio jets are typically weak, with Mach numbers <1.5 (For- 
|man et al^||2005| |McNamara et al.||2005| e.g.,). Shocks driven by 
the tips of the radio jets would have continued to move out into the 
ICM after the jets shut down, becoming increasingly difficult to de- 
tect as they move into regions of lower X-ray surface brightness. 
The northwest filament might represent a shock front, or merely 
material pushed out subsonically by an expanding cavity. Unfortu- 
nately the data are too shallow to allow us to determine whether the 
temperature of the northwest filament differs from the surrounding 
ICM. The strength of any potential shock can be estimated from the 
Rankine-Hugoniot jump conditions (e.g., |Landau & Lif shitz| 1 959) >, 
based on the density jump. Assuming temperature differences have 
a minimal impact on surface brightness, the density jump can be ap- 
proximated as the square root of the surface brightness jump. The 
surface brightness model does a rather poor job of representing the 
ICM in this region, but correcting for the underlying trend in sur- 
face brightness according to the model prediction, the remaining 
surface brightness jump is only a factor of 1 .08, suggesting a shock 
of Mach ~ 1 .03. This suggests that the northwest filament is likely 
to be the product of near-sonic expansion, rather than a supersonic 
shock. 

Estimating the age of the VLBA double is difficult, since its 
angular scale is far smaller than that of the X-ray pressure pro- 
file. We can calculate an approximate sonic timescale under the as- 
sumption that the central temperature (and therefore sound speed) 
is similar to that observed in our innermost projected temperature 
bin, ~3.3 keV. This implies a sound speed of ~870 km s _1 and 
an timescale of order ~10 J yr. The buoyant velocity depends on 
the gravitational mass inside the mean radius of the lobes, which 
is unknown, but would exceed the sound speed unless the enclosed 
mass is <2.4x 10 9 Mq, and would be ~2/3 of the sound speed for 
10 9 Mq. We have no way of recognising supersonic expansion on 
these scales, so the cavities could be younger than the sonic esti- 
mate. Conversely, we cannot confirm the presence of cavities, and 
if the lobes are partially filled by entrained thermal plasma the ef- 
fects of buoyant forces could be reduced to the point where the 
lobes remain in the core for long periods (as in NGC 5044, David 
|et al.|20 09). These estimates of the age of the VLBA lobes should 
therefore be treated with caution. 

As well as the age of the radio source, we are also interested 
in placing limits on the age of the current period of quasar activity. 
Current simulations of AGN in groups and clusters often model the 



feedback as either a bipolar outflow or a simple injection of thermal 
energy (e.g., |Dubois et al.|2010||McCarthy et al.|2010|[20TT) . Both 
cases produce larger-scale hot outflows which drive material into 
the outskirts of the system. The current radiative luminosity of the 
QSO (a few 10 47 erg s -1 ) is comparable to the mechanical power 
of the outflows modelled in these simulations, and the authors argue 
that radiative feedback should be efficiently converted into a me- 
chanical outflow in the group or cluster core ( |Dubois et al.|2010fr . 
If the QSO were producing a large-scale wind, we would expect it 
to inflate a central cavity and/or drive shocks out into the surround- 
ing ICM. We might expect to detect such a feature once it grew to 
a sufficient size to be resolved. Similarly, if the QSO were directly 
heating the ICM, we might expect to see a central temperature rise. 
We are limited by the bright central X-ray source, but we would 
expect to detect features of comparable size to our innermost tem- 
perature bin, 13.8 kpc (2.5")- For a sound speed of ~870 km s -1 
this suggests a timescale of ~ 15.6 Myr. The lack of features on this 
scale suggests that the AGN is either ineffective in heating the ICM 
or has only been powerful enough to do so for a relatively short 
time. It should also be noted that the QSO has as yet only heated 
the dust enshrouding it, rather than destroying or sweeping it out of 
the galaxy. 



7.4 The northeast optical filament and the radio spur 

The small residual X-ray feature northeast of the core in Fig. [7] 
extends in roughly the same direction as a spur of radio emission. 
Fig.[8]shows the 1.28 GHz radio and X-ray residual contours over- 
laid on the HST 622W-band image. The X-ray and radio features 
appear to be correlated with the brightest of the gaseous filaments 
north of the BCG. 

Optically luminous filaments of ionised gas have now been ob- 
served in numerous nearby cool core clusters, and they are known 
to be associated with the coolest, densest regions of the X-ray emit- 
ting ICM (e.g.,|Heckman et al.|1989||Conselice et al.|200Tj[Fabian 



et al.|2 003 , Cra wford et al . 2005 , Ha tch et al.||2007| |Hicks etaT 
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201 \\ . This association is clear in CL09, with the optical filaments 
falling within the coolest part of the temperature map. This sug- 
gests that the ionised gas is either the product of cooling from the 
ICM, or that it enhances cooling through conduction of heat from 
the hot gas or by mixing with the hot phase. 

A correlation with radio emission is more difficult to explain. 
The spur is detected at ~6<r significance in the 1.28 GHz map, 
and a similar feature on the far side of the central component is 
detected at ~3<r significance. The northeast spur should be detected 
in the VLA 1.5 GHz map, but is not; however a feature coincident 
with the southwestern counterpart is detected. We cannot be sure 
whether these spurs are real structures. A number of other small 
structures are seen in the 1.28 GHz map which seem likely to be 
noise features or marginally detected background sources (e.g., the 
small blob of emission attached to the tip of the northern lobe). 
We therefore consider their possible origin, but note that deeper 
observations are needed to confirm or refute their existence. 

Star formation is observed in some optical filaments in galaxy 
clusters (e.g.,|Johnstone et al.|1987||McNamara & Q'Connell|1993[ 
Voit & Donanue||19971|Crawford et al.||19991 |Q'Dea et al.||2004l 
Canning et aTl[20T0l |Oonk et al.||20TTj |McDonald et^[2fjTl) , 
and could produce some degree of radio emission. Measuring the 
1.28 GHz flux density of the spur is difficult, given its proximity to 
the bright central component, but our best estimate is ~0.25 mjy. 
Using the known relation between radio flux and star formation 
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rate jBeiIl|2003l |Garn et al.||2009) , we estimate that an SFR of 
~95 Mq yr~ would be required to produce the feature. Star for- 
mation at this level would likely be obvious in the optical images 
and wou ld have been detecte d in previous spectroscopy. For com- 
parison, l Bildfell et al. 2008 ) estimate an SFR or ~41 M e yr~ 1 for 
the whole BCG. 

The spurs could be evidence of a third epoch of jet activity, 
with a different axis to both the large-scale jets and the VLBA 
double. If the spurs represent the full extent of a jet, we might ex- 
pect small-scale cavities at the position of the spurs. It is unclear 
whether we could detect such small features, but the bright X-ray 
feature associated with the northeast spur would certainly not be 
expected. 

Alternatively, if the spurs were remnants of a jet which was 
active before the large-scale jets but on a similar spatial scale, their 
radio emission at larger radii might have faded beyond our abil- 
ity to detect. In this case the optical filaments might be gas up- 
lifted from the core by cavities which have now moved out to large 
radii where they cannot be detected in the Chandra observation. 
A similar correlation is seen in the Perseus cluster, between opti- 
cal filaments behind an old cavity and a ridge of emission in the 
radio mini halol jFabian et al.|201 1) , but we note that there is no ev- 
idence of any diffuse radio halo in CL09 down to the sensitivity of 
our low-frequency observations. The suggestion that the spurs are 
a remnant of a previous outburst is thus highly speculative and the 
data provide no supporting evidence for this hypothesis beyond the 
presence of the spurs. 

Only a deep higher resolution observation can determine 
whether the spurs are real and physically correlated with the other 
structures in the cluster core. We therefore conclude that while the 
X-ray/optical correlation is related to cooling in the dense ICM, 
the reality and origin of the radio emission is uncertain and cannot 
be resolved with the current data. 



/ 

/ Radio 



Ionization 



X-ray 



Galaxies 



4" / 22 kpc 



Inner Double 



Figure 18. HST 622W filter image of the core of the core of CL09. The 
axes of the old radio jets (red dashed line), inner double radio source (blue 
solid line), major axis of the cluster X-ray emission (from our fits to the 
Chandra data, solid black lines), major axis of the cluster galaxy distribution 
(black dashed lines, Hines & Wills 1 993 1 and the apparent opening angles 
of the AGN torus (solid black arcs), estimated from the dilution-corrected 
polarized optical flux (Hines et al.| 1 99 9 1, are indicated schematically. The 
position (but not the angle) of the radio axes has been adjusted so that they 
pass through the peak of optical surface brightness, to allow direct compar- 
ison with the torus opening angles. The offsets are <0.5". 



8 DISCUSSION 

8.1 History of the AGN and the cluster 

Previous studies of CL09 have noted several features which suggest 
that the orientation of the AGN, radio jets and surrounding dusty 
torus may have recently changed. These include the misalignment 
of the radio jets and the ionisation cone of the quasar (as traced in 
polarized optical line emission, Hines & Wills 1993, Crawford & 
|Vanderriest|T996l|Hines et al.|1999||Armus et al.|1999||Tran etaT*| 
|2000[ >, the steep high-frequency spectral index of the radio lobes, 
and the flatter index and high luminosity of the radio core found 
by |Hines & Wills ( 1993}. It was argued that the large-scale radio 
jets were no longer powered by the AGN, and that some event had 
caused the AGN to change its axis and launch new jets at an angle 
close to the line of sight, which were as yet too small to be resolved 
( |Hines & Wills|1993||Crawford & Vanderriest|1996) . 

Our results support this general picture and improve on it. 
Fig. [T8] shows in diagrammatic form the orientation of the various 
structures associated with the AGN and BCG. As noted by |Hines| 
|& Wills | H993[ l the axis of the large-scale radio jets falls outside 
the ionisation cones of the AGN. The VLBA inner double source 
has a position angle ~80° (from west), within the opening angle 
of the ionisation cones and torus. However, while it was previously 
argued that any new jet must be inclined close to the line of sight, 
the similar sizes of the two components of the VLBA inner double 
suggest that it is in fact aligned close to the plane of the sky. While 



the projected orientation of the VLBA source is consistent with the 
ionisation cones, they do not appear to be perfectly aligned, and an 
orientation close to the plane of the sky would conflict with the in- 
clination of ~37°to the line of sight found from optical polarisation 
measurements ( |Hines etal.|i999"l >. This may indicate that as well as 
a change in AGN orientation between the epoch of the large-scale 
jets and that of the VLBA double, there has been a further change 
and the inner double has shut down. 

The current state of the radio jets is still open to question. 
Both dynamical and synchrotron ageing arguments suggest that the 
large-scale jets are old, and their spectral index confirms that they 
have not been powered by the AGN for some considerable time. 
The steep spectral index of the central component of the large- 
scale source suggests that the inner double is also old. In this case 
the lack of a resolved core in the VLBA data and the misalignment 
of the inner double with respect to the ionisation cones would be 
explained; the AGN would be in a radiatively efficient accretion 
state with no current radio activity. However, given the small size 
of the inner double (~200 pc) and short dynamical timescale for 
the lobes (~10 5 yr), and the lack of a direct measurement of their 
spectral index, it is also possible that the core is active but self- 
absorbed at 1.4 GHz. Multi-frequency VLBA observations would 
be needed to measure the spectral index of the inner double and 
resolve this question. 

|Hines & WilIs| ( |T993| > also noted that the radio axis is roughly 
aligned with the minor axis of the distribution of galaxies within 
CL09 and suggested that this correlation could be driven by the 
pressure gradient of the ICM. If buoyant forces are important in 
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the growth of the radio jets, they will tend to align themselves with 
the steepest pressure gradient, in this case along the minor axis of 
the cluster ICM. Fig. [T5] shows that both the large-scale radio jet 
and the VLBA inner double are aligned close to the minor axis of 
the cluster galaxy distribution, but are somewhat less well aligned 
with the ICM. The alignment may be coincidental, as the pressure 
gradients on 200 pc scales are probably trivial and therefore the 
alignment of the VLBA double is unlikely to be governed by buoy- 
ant forces. It should also be noted that the coolest part of the ICM 
is more clearly perpendicular to the jet axis, but this is probably 
caused by the excavation of cavities by the jets, rather than the tem- 
perature distribution affecting jet propagation. 

More interestingly, the various AGN axes are poorly aligned 
with the BCG, whose major axis runs almost north-south. |Craw-| 
|ford & Vanderriest] fT?9 6j) note that the velocity structure of the 
[Om] emission in the BCG suggests rotation about a roughly east- 
west axis. The relatively small rotational velocity (<200 km s _1 ) 
argues against the gas being located in a simple equilibrium rotat- 
ing disk given the deep potential well of the BCG, but a rotating 
structure which has yet to reach equilibrium with the galaxy seems 
possible. The coincidence between the onset of the most recent pe- 
riod of star formation (<200 Myr ago Pipino et al. 2009 ) and the 
cessation of jet activity (~130 Myr ago), both of which could be 
caused by an increase in the available cool gas supply, suggests 
that there was a sudden influx of gas into the galaxy ~200 Myr ago. 
This raises the possibility that the BCG has recently undergone a 
gas-rich merger. The presence of dust mixed with the ionised gas 
in the inner filament and galaxy core ( [Tran et al.|2000) supports 
this hypothesis. Such a merger could provide the cold gas needed 
to fuel the current powerful star formation and rapidly increase the 
accretion rate of the AGN, causing it to shift from a radiatively in- 
efficient mode with large-scale jets, to a radiatively efficient quasar 
mode. The axis of the AGN could also have changed in response 
to the different angular momentum of the infalling material, and 
might not yet have stabilised. 

CL09 is a massive cluster with an extended dense ICM. A sin- 
gle gas-rich galaxy falling into the cluster would likely be stripped 
before it could merge with the BCG. However, a gas-rich galaxy in 
the core of an infalling galaxy group or cluster could be protected 
by the surrounding intra-group medium, since this would need to be 
stripped before galactic gas could begin to interact with the ICM of 
CL09. This could potentially insulate the galaxy until the two clus- 
ter cores merged. Simulations suggest that a cluster-cluster merger 
beginning ~4-5 Gyr ago could produce an apparently relaxed clus- 
ter with weak spiral disturbances in the ICM, but with the clus- 
ter cores merging only in the last few hundred Myr (Poo le et al.| 
|2006[ >. A merger scenario therefore seems plausible. The distribu- 
tion of galaxies could potentially provide evidence of a merger, and 
it is intriguing to note that an apparent chain of galaxies extends 
from northeast of the BCG across its northern side, overlapping 
the northern optical filament. However, without velocity measure- 
ments from these galaxies, we cannot be certain that this is not a 
mere chance superposition. 

Assuming that the star formation is taking place in a disk and 
follows the standard Kennicutt-Schmidt relation (Kennicutt 1998), 

Af* = 0.017 x M g asVc/R (2) 

we can use the star formation rate (M*=41 Mq yr~ 
|fell et al.||2008) and available gas mass ( |Combes et al.||2011 



M a 



-^3.2x 10 Mq) to estimate the radius of the star formation 



region, R. Takin g the rotational velocity of the [Om] emission 
(u c ~200 km s" 1 [Crawford & Vanderriestl|l996| to be represen- 



tative, we find that the disk is likely to have a radius _R~270 pc 
(~0.05"). This radius is considerably smaller than the region of 
blue optical colours in the core of the BCG (~20 kpc radius, |BlkTl 
|fell et al.|2008} . The blue core is probably produced by a relatively 
small mass of young stars formed during a more energetic period of 
star formation shortly after the cold gas entered the BCG ~ 100 Myr 
ago. The ongoing star formation traced by the [Oil] emission is 
probably fuelled by material left over from this initial burst, which 
has had time to fall into the central regions of the BCG. A com- 
pact, dense, gas-rich star forming disk around the AGN provides 
a natural explanation for the heavy obscuration seen in X-ray and 
IR-NUV SED fitting (e.g., |Vignali et al.|2011| l. 

The presence of cavities, and in particular their correlation 
with the jets but not the lobes, complicates our picture of nuclear 
activity. For the southern jet it seems possible that the observed 
cavity is part of a larger structure extending out to include the lobe, 
but that the Chandra observation is not deep enough to detect its 
outer parts. However, if the northern cavity was formed by the jet, 
the rim of material around the cavity must have been formed at the 
same time, with the northern jet first inflating the cavity and then 
pushing through its outer rim to form a lobe outside. Again, the 
northern lobe would in this case have an associated cavity, which 
we fail to detect because the observation is too short and the lobe 
is partly obscured by the rim of the inner cavity. This formation 
of multiple cavities might be explained by variation of jet power 
during the outburst, or by pre-existing structures in the ICM. 

An alternative scenario is that the cavities were formed by 
a separate AGN outburst, and are unrelated to the currently visi- 
ble radio emission. If the cavities were formed first, by some ear- 
lier AGN outburst, the weak correlation between radio and X-ray 
structure would be explained, since the jets need not pass through 
the cavities. We would require only that the jet axis changes with 
time piabul et al.||2012"| and references therein), or that motions 
within the ICM have moved the cavities away from their original 
axis (IMorsony et al. 2010). The opposite case, in which the jets 
formed first and the cavities through a second outburst, would re- 
quire the radio emission from this second jet to be disguised. It is 
just possible that the apparent knots in the radio jet are in fact the 
lobes of some later outburst. However, this seems unlikely given 
their poor correlation with the cavity morphology and the consis- 
tent spectral indices across the jets and lobes. 

In either case, the straightness of the radio jets and the fact 
that they extend across the boundaries of the spiral X-ray residu- 
als strongly suggests that the position of the BCG relative to the 
surrounding ICM has remained fairly constant over the last ~100- 
150 Myr. If the BCG position has remained relatively static, any 
cluster-cluster merger must have occurred several Gyr ago, al- 
lowing time for the cluster to have largely relaxed. It also seems 
likely that any infalling system which provided the gas to fuel star- 
formation and the AGN was significantly less massive than the 
BCG. 



8.2 The QSO as a source of feedback heating 

Energetically, it seems clear that the radio jets are capable of bal- 
ancing cooling in CL09, provided that: i) the duty cycle is relatively 
long, so that there are not long periods of cooling with no jet ac- 
tivity, and ii) that the energy injected by the jets can be efficiently 
coupled into the ICM. These are conditions which apply to all sys- 
tems in which feedback by AGN jets has been proposed, though 
in CL09 the evidence of multiple epochs of jet activity make it 
more likely that condition (i) is being met, at least over timescales 
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~50-150 Myr. The presence of cool gas and star formation in clus- 
ter cores is thought to suggest that cooling from the hot ICM is 
not been completely stopped by AGN feedback (e.g., Bildf ell et al.| 
|2008| |Pipino et al.|2009) . In CL09 we only require cooling to be 
suppressed and, since the mechanical energy involved in inflating 
the cavities (PV) is roughly equal to the cooling luminosity of the 
cluster core, it seems that the jets are easily capable of providing 
this suppression. It should also be noted that gas may have been 
brought into the cluster core via a merger, rather than cooling from 
the ICM. 

There is no evidence of heating of the hot ICM through ra- 
diative feedback. T he accretion related luminosity of the AGN (a 
few x 10 47 erg s~ 1 , |Vignali et al.|201 l|[Ruiz et al.|20"To] > is a fac- 



tor ~100 higher than our estimated jet powers, and comparable to 
the Eddington luminosity for a 1.5 x 10 9 Mq black hole, similar to 
those in nearby BCGs such as M87, or ~10 per cent, of the Ed- 
dington luminosity if the black hole mass is of order 10 10 M©. It 
therefore seems likely that the AGN is close to its Eddington lumi- 
nosity, and even with a lower coupling efficiency than for jet feed- 
back we might expect to see some impact on the ICM. Radiative 
heating of the ICM might occur directly, or via a wind or radiation- 
driven convection. Optical spectra of the AGN show a blueshifted 
component with a velocity offset of ~1300 km s _1 (Crawford & 



Vanderriest 1996) suggesting an outflow at least on small scales 
However, if the ICM is being heated we might expect to see either 
a central temperature rise, a central cavity, or structures associated 
with outflowing gas. There is no evidence of any such structures, 
at least on the 2.5" (~14 kpc) scales we can resolve outside the 
central point source, and thus no evidence of heating on the scale 
needed to balance cooling. 

This lack of evidence could be an issue of power or of 
timescale. It has been suggested that photoionization and Comp- 
ton heating may only be able to raise the temperature of gas to 
~1 keV {Sazonov et al.|2003] i. If radiative feedback can heat gas 
above this temperature, following Siemiginowska et al. (2010) we 
can estimate the efficiency of the coupling between the QSO radia- 
tion and the surrounding ICM. For the cool core, the opacity of the 
gas to Compton scattering is r = n e roaT where n e ~ 0.05 cm -3 
within a radius ro=70 kpc, and <tt is the Thompson scattering 
cross-section. This suggests that ~0.7% of photons will interact 
with gas in the core, but only the fraction of these with energies 
>1 keV could heat the gas. This fraction is not well known, but 
is likely to be <10%. Taking into account these factors, the lumi- 
nosity available to heat the gas may only be a few x 10 44 erg s~ , 
comparable to the cooling luminosity. The quasar could thus be too 
weak to stop cooling, despite its apparent high luminosity. Alterna- 
tively, the QSO may simply not have been in its current state long 
enough to affect the cluster core. If our hypothesis of a gas rich 
merger is correct, the transition to a quasar state was triggered by 
an increase in accretion rate which also shut down the radio jets, 
and is fuelled by the same material which fuelled star-formation. 
Observational studies of galaxies hosting moderate-luminosity su- 
permassive black holes have found a delay of ~ 100 Myr between 
the times of greatest star formation activity and peak AGN lumi- 
nosity (Schawinski et al. 2009). This is supported by simulations, 
which suggest that the delay is caused by the need to remove angu- 
lar momentum from cold gas brought into the galaxy via mergers 
(Hopki nspOT2) . IRAS 09104+4109 is a considerably more massive 
system than those considered by these simulations, but a similar de- 
lay is certainly possible. The QSO has had at most 100-200 Myr to 
heat its surroundings, and a dynamical delay in fuelling could mean 
that it may only have reached its current high radiative power more 



recently. However, given the lack of any clear effect on the ICM, we 
must conclude that in CL09, there is only evidence of feedback via 
the radio jets, and that only this feedback mode is required to bal- 
ance cooling from the hot ICM. In this regard IRAS 09104+4109 
appears similar to HI 82 1+643 ( |Russell et aLpOTo) and 3C 186 
(Siemiginowska et al. 2010), both (unobscured) quasars located in 
cooling flow clusters whose radiative output appears to have had 
little impact on their environment. 



9 CONCLUSIONS 

Using a combination of deep X-ray and multi-frequency 
radio observations, we have examined the type-II QSO 
IRAS 09104+4109, located in the BCG of the rich, cooling 
flow cluster CL 09104+4109. Prior studies have shown that the 
BCG is an HLIRG powered by a combination of AGN emission 
and star formation, and hosts an old, steep spectrum FR I-FR II 
radio source which has inflated cavities in the ICM. Our analysis 
allows a more detailed examination of the relationships between 
the cluster, BCG and AGN, and provides insight into the history of 
the system. Our results can be summarised as follows: 

(i) We present new GMRT observations at 240, 317, 610 and 
1280 MHz, and reanalyse archival VLA observations at 1.5 and 
4.8 GHz. We confirm the previous findings that the radio source 
has symmetrical straight jets ending in small lobes, with a bright 
central component. The radio spectrum of the source is an unbro- 
ken powerlaw extending down to <240 MHz, with a steep spectral 
index a=1.25±0.01. Spectral index variations across the source are 
minimal, and the central source has a comparable index to both the 
northern and (within the uncertainties) southern lobes. The steep 
spectral index indicates that the jets are not currently active and the 
source is passively ageing. The straightness of the jets suggests that 
the position of the BCG relative to the local ICM has not signifi- 
cantly changed over the lifetime of the source. 

(ii) Analysis of a previously unpublished archival 1.4 GHz 
VLBA observation of the central source resolves it into an extended 
double source ~200 pc across. Both components are extended, 
suggesting that these are small-scale inner radio lobes, probably 
formed by a separate period of jet activity to the large-scale jets. 
The dynamical timescale of the lobes is ~10 5 yr. The inner dou- 
ble may therefore be young compared to the large-scale radio jets, 
unless the buoyancy of the lobes is suppressed by mixing with en- 
trained thermal plasma. However, while the spectral index of the in- 
ner double cannot be measured directly, its flux is identical (within 
uncertainties) to the flux measured for the central component of the 
large-scale source and there is no evidence of significant additional 
emission which might affect the spectral index measured on larger 
scales. If the spectral index of the central component is an accurate 
measure of the index of the inner double, the inner lobes are not 
currently powered by the AGN. No radio core is detected, but it is 
unclear whether this is because the AGN is not currently launching 
jets, or whether the core is active but self-absorbed at 1.4 GHz. The 
axis of the inner double differs from that of the large-scale jets by 
~17°, supporting the hypothesis the AGN axis has changed. 

(iii) We measure deprojected radial profiles of temperature and 
density of the ICM out to ~850 kpc, and determine the gravita- 
tional mass and gas mass profiles of the cluster. The core temper- 
ature of ~4 keV and peak temperature of ~8 keV agree well with 
previous estimates ( Hlavacek-Larrond o et al.|2011) . The total mass 
profile agrees well with gravitational lesing measurements, and the 
cluster falls on both the L-T and M-T relations. The central cooling 
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time within 5" (~27.5 kpc) is 9.9±0.9xl0 8 yr, and we estimate 
an isobaric cooling rate of <235 Mq yr -1 within this radius. We 
map the 2-dimensional temperature structure of the cluster core, 
and find that the coolest gas (kT<5 keV) has an asymmetrical dis- 
tribution anti-correlated with the axis of the large-scale jets. The 
cool, optically-luminous gas filaments detected by HST lie within 
this region of cool ICM temperatures. 

(iv) We confirm the presence of X-ray cavities northwest and 
southeast of the BCG, coincident with the radio jets, at a radius 
~24 kpc. However, the precise relationship between the jets and 
cavities is unclear. The northern cavity is bounded on its west and 
northwest side by regions of excess surface brightness which may 
represent a shell of gas expelled from the cavity. The southern cav- 
ity appears larger, extending along the jet and to the east outside 
the radio contours, and lacks any clear shell. Cavities are not de- 
tected coincident with the radio lobes, but the Chandra exposure is 
probably insufficient to detect cavities of the expected volume. We 
estimate the enthalpy of the detected cavities (~7.7x 10 60 erg for 
both, summed) and for any undetected cavities at the location of 
the radio lobes. The northwest rim of the northern cavity has a rel- 
atively sharp, linear appearance, raising the possibility that it may 
be a shock. Unfortunately the Chandra observation lacks the depth 
to detect any temperature increase above that of the surrounding 
ICM, but based on the surface brightness increment we find that if 
the feature is a shock it is only mildly trans-sonic. 

(v) We find that the large-scale jets are no longer powered 
by the AGN, but are relics of activity which ceased ~ 150 Myr 
ago. The radiative age estimate derived from the radio spectral in- 
dex {t r ad ~130 Myr) is in agreement with the refill timescales 
estimated for the cavities associated with the radio jets (~120- 
160 Myr). The radiative and dynamical timescales are in agreement 
with optical/NUV estimates of the likely age of the star formation 
in the BCG, ~70-200 Myr. The similarity between these timescales 
and the presence of dust in theionised gas of the inner filament, and 
in the galaxy core, suggests that the galaxy experienced a sudden 
influx of cool gas ~200 Myr ago, most likely as a result of a merger 
with a gas-rich interloper. A rapid inflow of gas from such a merger 
could have triggered star formation and increased the accretion rate 
of the AGN, causing it to cease jet production and shift into a ra- 
diatively efficient accretion mode. 

(vi) In both Chandra and XMM data, subtraction of the best 
fitting elliptical surface brightness /3-model leaves a spiral resid- 
ual feature in the ICM, extending to ~200 kpc. This is similar to 
features seen in other clusters and in simulations, associated with 
"sloshing" motion of the cluster core caused by a tidal encounter 
with another massive system, or with a previous merger. No sharp 
surface brightness or temperature features are observed, and the 
temperature map shows no correlation with the spiral residual, sug- 
gesting that the cluster has had time to relax. A merger with a 
galaxy group or cluster, beginning ~4-6 Gyr ago with core merger 
occurring 0.5-1 Gyr ago, might explain both this spiral feature and 
the evidence of a recent gas-rich galaxy merger with the BCG. 

(vii) We estimate the mechanical power output of the radio jets 
from the cavity enthalpies. Assuming that the jets are efficiently 
coupled with the ICM, we find that the jets are capable of balancing 
radiative cooling within the radius to which they extend (77 kpc), 
and probably out to at least 200 kpc. Cooling times at these radii are 
>9.5 Gyr, and it therefore seems likely that the jets are capable of 
maintaining the thermal balance of the cluster as a whole. While no 
direct evidence of gas heating by the AGN jets is observed, weak 
shocks of the type usually associated with jet sources would now 



be difficult to detect, given the age of the outburst and the relatively 
shallow Chandra data. 

(viii) Although optical spectra of the AGN may suggest an out- 
flow on small scales, we see no evidence of a large-scale radiatively 
driven wind of the sort predicted by some simulations of quasars in 
dense environments. The lack of any central cavity, temperature 
rise, or evidence of outflowing gas suggests that either the quasar is 
insufficiently luminous, that it has had insufficient time to produce 
detectable changes, or that radiative feedback is ineffective in heat- 
ing the hot ICM. It is also notable that the quasar is still enshrouded 
by ~3 x 10 9 Mq of cool gas (inferred from CO measurements) and 
<2xl0 7 Mq of dust. While some of this material has been heated 
and/or ionised, the quasar has not destroyed it or swept it out of the 
BCG. The estimated bolometric luminosity of the quasar is a few 
x 10 47 erg s _1 and if our merger hypothesis is correct it entered its 
current radiatively efficient state ~50-150 Myr ago. 
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